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ABSTRACT 

We investigate the observational requirements for the detection of sound-wave-like features in 
galaxy cluster cores. We calculate the effect of projection on the observed wave amplitude, and 
find that the projection factor depends only weakly on the underlying cluster properties but 
strongly on the wavelength of the sound waves, with the observed amplitude being reduced by 
a factor ~ 5 for 5 kpc waves but only by a factor ~ 2 for 25 kpc waves. We go on to estimate 
the time needed to detect ripples similar to those previously detected in the Perseus cluster in 
other clusters. We find that the detection time scales most strongly with the flux of the cluster 
and the amplitude of the ripples. By connecting the ripple amplitude to the heating power in 
the system, we estimate detection times for a selection of local clusters and find that several 
may have ripples detected with ~ 1 Ms Chandra time. 
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1 INTRODUCTION 

The thermal history of the hot gas in galaxy clusters has been a 
long-standing puzzle to astronomers. In many clusters (often re- 
ferred to as cool-core clusters), the cooling time of the hot X-ray 
emitting IntraCluster Medium (ICM) drops to a fraction of the age 
of the universe in the brightest central regions. The first models of 
these clusters suggested that a significant amount of the gas should 
form a "cooling flow", in which gas in the cluster core cools out of 
the X-ray band, reducing the pressure support on the surrounding 
material and causing the outer gas to flow toward the centre. The 
cool gas predicted by this model is expected to fuel star formation 
in the brightest cluster galaxy. However, it has long been known 
that the star formation rates in these systems are significantly lower 
than these cooling flow models would predict. 

With the launch of the most recent generation of X-ray 
satellites — particularly Chandra and XMM-Newton — it has 
been shown that the amount of cool X-ray emitting gas in clus- 
ters is much smaller than p redicted by the simple cooling flow 
models jPeterson et all [2003). This has led to the widely ac- 
cepted conclusion that the ICM is being heated, reducing the 
coo ling flow to around 10 p ercent of the expected valu e (see 
e.g. lPeterson & Fabianll2006l or lMcNamara & NulseJl2007l for re- 
views). 

At the same time that Chandra and XMM-Newton have con- 
firmed that heating of cluster cores is required, there have been 
substantial advances in identifying possible heating mechanisms. 
Whilst there are a large number of possible heating mechanisms, 
the current leading contender is mechanical energy injection by the 
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central black hole. Evidence of the interaction between the ICM 
and the central clu ster was first noted in ROSAT observations of 
the Perseus cluster (B ohrin g er et al . ; 1993), with relativistic plasma 
from the AGN jets displacing the ICM to form radio-emitting cav- 
ities. Since that time, detection of such cavities in nearby clusters 
has be come commonplace, with a recent study bv lDunn & Fabiar] 
(2006) showing evidence for X-ray cavities in at least 70 percent 
of the clusters in their sample with short central cooling times and 
central temperature drops. The work done on the ICM in the in- 
flation of these cavities provides a plausible source for the heating 
of the ICM, and several studies have shown that the power associ- 
ated with the bubbles' inflation is sufficient to offset the cooling lu- 
minosity jBirzan et al J 120041 : iDunn & Fab"iai]|2004 iRaffertv et all 
l2006l : lDunn e"tai]|2005l) . 

Despite the promise shown by the AGN-injection model of 
cluster-core heating, there are a substantial number of details that 
remain unclear. One particular problem is how the energy injected 
by the AGN near the cluster core is distributed over the whole cool- 
ing region in a quasi-isotropic manner. The peaked central metal 
abundances in many cool-core clusters indicate that this process 
must be gentle, suggesting that it is not simply a matter of clus- 
ters having occasional unusually large outbursts that heat the outer 
regions. 

A possible clue to the method of energy transport away from 
the cluster core was given by the unexpected detection of quasi- 
spheric al ripples in the sur face brightness structure in the Perseus 
cluster jFabian et alj|2003l) , shown in Fig. [TJ with a high pass fil- 
ter app lied to remove the un derlying cluster emission (reproduced 
from Sanders & Fabianl2007l) . Ripples in surface brightness are in- 
dicative of fluctuations in the density of the ICM and consequen- 
tially, these features have been interpreted as sound waves gener- 
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Figure 1. A 900 ks Chandra image of the Perseus cluster core, with a 
high pass filter applied to remove the underlying surface brightness pro- 
file and show the ripples. It is clear that the ripples are highly symmetric 
about t he central source, but not coherent over the full 2tt. Reproduced 
from lSanders & Fabianl J2007I) 
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Figure 2. Amplitude of surface brightness ripple height as a percentage 
of the density ripple amplitude for a range of underlying density model /3 
parameters and density perturbation wavelengths and with fixed amplitude 
density perturbation / = 0.15 and fixed density model r core = 25 kpc. 



ated by the cavity inflation process. Such sound waves are expected 
to be effective in distributi ng the AGN's mechanical e nergy input 
over the entire cluster core 1 Mc Namara & N ulsen 2007) and, in the 
case of Perseus, it has been shown that the en ergy they carry is 
compa rable to that required to offset the cooling lSanders & Fabianl 
( 120071) . 

Whilst cavities are common, associated sound waves have 
never been detected in a cluster other than Perseus. However, the 
Perseus cluster is both a factor of ~ 2 brighter than any other clus- 
ter and has a factor 2 more observation time than any other clus- 
ter. In order to understand whether sound waves do indeed play 
a general role in the heating of cluster cores, however, their pres- 
ence needs to be confirmed in other systems. Our main aim in this 
paper is to investigate the relationship between the intrinsic and 
observed properties of sound waves in galaxy clusters and to as- 
sess the prospects for detecting sound waves in systems other than 
Perseus. To do this, we will not try to model the detailed time evo- 
lution of ripples in cluster cores, but will assume that the almost- 
monochromatic, low amplitude, ripples in Perseus are a good tem- 
plate for the features expected in other clusters. By imposing rip- 
ples derived from this template on model clusters, we will deter- 
mine how the detection time varies with the properties of the ripples 
and the underlying cluster profile. 

Throughout the paper we adopt a cosmology where Hq = 
70kms~ 1 kpc -1 , Om = 0.3 and Qa = 0.7. The redshift of the 
Perseus cluster is 0.018. 



2 EFFECT OF PROJECTION ON RIPPLE AMPLITUDES 

Inferring information about extended sources such as galaxy clus- 
ters is complicated by the fact that, whilst the source is intrinsi- 
cally three dimensional, our view of it is only two dimensional. In 
the case of an idealised spherically symmetric cluster, the image as 



viewed on the sky is related to the underlying emission profile by 
the Abel integral: 
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where b is the projected distance on the sky, r is the 3-dimensional 
radius, e is emissivity and I is the surface brightness. 

Ignoring the effect of projection, a monochromatic density 
wave in an isothermal galaxy cluster with e a will lead to an 
emissivity profile of the form 
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where / is the fractional amplitude of the density wave. The corre- 
sponding amplitude of the emissivity fluctuations is then 
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i.e. the amplitude of the rip ples in emissivity is \h s \ ~ 2/. 

However, as noted by 'Fabian et al.l d2006t) . the effect of pro- 
jection is to reduce the magnitude of the surface brightness ripples 
so that in general \h s b\ < \h e \, with troughs in the emissivity being 
filled by the surroundi ng peaks, an d vice- versa. For the case of the 
Perseus cluster, dFabian et all2006T) found that projection acts such 
that \h s b\ ~ //2.5 ~ |/i e |/5. However this result is not universal 
as the extent to which a trough in emissivity is filled by the sur- 
rounding peaks depends on the relative overlap and emission of the 
peaks and troughs along the line of sight. This is a function both of 
the underlying surface brightness profile and the properties of the 
ripples themselves. 

To quantify the relationship between the wavelength and clus- 
ter properties and the projected ripple height, we have numeri- 
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Figure 3. Amplitude of the surface brightness ripples as a percentage of 
the density ripple amplitude for a range of ripple wavelengths and beta- 
model core radii for fixed amplitude density perturbation / = 0.15 and 
fixed p = 0.6. The amplitude is almost independent of core radius for 
wavelengths below 25 kpc. 



cally calculated the projected ripple height for monochromatic rip- 
ples in clusters with an underlyin g beta-model density distribution 
dCavaliere & Fusco-Fe miano 197a): 
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where r c is the core radius of the cluster, f3 is the beta parameter 
and n e fi is a normalisation. The beta model provides a reasonable 
fit to the density distribution in many clusters except in the inner- 
most region where adding a second beta-model component often 
provides a better fit. 

Fig|2]shows the median projected ripple height calculated for 
a isothermal kT = 6 keV cluster with an underlying beta-model 
density distribution, 
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with values of (5 in the range 0.2 — 1.8 and monochromatic den- 
sity perturbations with / = 0.15 and wavelengths in the range 
5 — 50 kpc. The use of the median height reflects the fact that not 
all ripples have exactly the same height in projection, however the 
heights are consistent enough that the exact method of quantifying 
|/i s i,| for a given set of parameters does not significantly affect the 
results. The projected height decreases with shorter ripple wave- 
lengths and shallower density profiles, which is expected as these 
are the conditions in which adjacent peaks and troughs in emissiv- 
ity will have the most similar amplitude along the line of sight. Fig. 
[3]shows the projected height for a similar cluster but with f3 = 0.6 
and a range of core-radii. It is apparent that the amplitude of the 
surface brightness ripples is not a strong function of the core ra- 
dius, although larger core radii — flatter emissivity profiles — do 
tend to show smaller surface brightness ripples, as expected. 

Testing clusters with fixed f3 and core radius, but varying den- 



sity amplitude, shows a linear variation in the projected amplitude 
over the entire range of density amplitudes likely to be relevant in 
clusters. This means that our calculations for a single ripple ampli- 
tude scale in a simple way to all amplitudes. 



3 TIME TO DETECT RIPPLES IN CLUSTERS 
3.1 A Naive Approach 

The basic condition to detect a ripple is that the variation in the 
number of counts in the peaks and troughs of the surface bright- 
ness profile are distinguishable from variations due purely to noise. 
Assuming a regime where Poisson noise dominates, this implies 
that to detect a single peak we require an observing time such that: 
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where n r i P pie is the average count rate in the ripple peak, n nor ippie 
is the average count rate that would be observed in the peak with- 
out a ripple, t is the time to detect the ripple, and a is the number 
of standard deviations above the mean we require for a detection. 
Therefore 
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Taking the underlying density profile to be a beta model and the 
cluster to be isothermal, the surf ace brightne ss profile will be a pro- 
jected beta model, of the form (Ettori 2000): 



S(r) = So (l + (r/rc) 



2x0.5-3/3 



(9) 



If the cluster is perfectly spherically symmetric and ripples 
are sinusoidal with a wavelength A and constant fractional surface 
brightness amplitude h the detection time is: 
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and S(r) — So&(r) represents the projected surface brightness in 
counts per second per unit area at radius r. ro, the lower limit of 
integration, is an integer number of wavelengths from the centre 
so the surface brightness is integrated over a peak. Taking some 
parameters appropriate for a Perseus-like cluster at z ~ 0.018; 
f3 = 0.5, r c = 30 kpc, A = 15 kpc, h = 0.05, So = 
5 counts s _1 arcmin -2 and a = 3, we find tdetect = 1.5 ks. This 
is clearly short compared to the time actually required to detect the 
sound waves in Perseus, which were not observed in a 25 ks obser- 
vation but were observed in a 200 ks observation. Part of this dis- 
crepancy can be explained by the fact that our calculation assumes 
that the ripples are coherent over a entire annulus of the cluster. In 
practice, the ripples in Perseus are coherent over a much smaller 
angle. To detect a ripple over some fraction \ of an angle requires 
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that the integration time is increased by a corresponding factor of 

i/x. 

Although this method is simple and easy to apply, the assump- 
tion of a uniform monochromatic wave may lead to misleadingly 
small detection times. For this reason, we have developed a method 
for estimating the detection time numerically for a more general 
waveform. 

3.2 An Improved Approach 

To better constrain the time required for detection of ripples in dis- 
tant clusters, we have developed a simple algorithmic approach to 
determining whether a cluster profile contains ripples. The input to 
the algorithm is a count rate profile either obtained from data or, for 
our purposes in the current work, generated artificially according 
to a method that will be described below. By generating artificial 
count rate profiles, we will be able to determine how the detection 
time for the ripples depends both on properties of the underlying 
cluster atmosphere and properties of the ripples themselves. 

The algorithm we use to determine whether a cluster has de- 
tectable ripples is: 

• Obtain a count rate profile R, either from data or from a model 

• Optionally (for simulated profiles) add Poisson noise to the 
count rate (assuming the input profile is an average count rate) 

• Convert the count rate profile to a surface brightness profile 
using the area of each annulus in the profile, to give an actual sur- 
face brightness profile S. 

• Fit a projected beta model to the surface brightness profile to 
give an underlying surface brightness profile B. 

• Calculate the fractional residuals of the actual surface bright- 
ness profile compared to the input surface brightness profile a = 
(S-B)/B 

• Fit a ID smoothing spline to the surface brightness residuals a 
to give a smooth profile a s . The degree of smoothing of the spline 
was adjusted by eye to provide an acceptable compromise between 
over-smoothing the ripples and closely following the noise in 5*. 

• Use the zero points of a s to delimit peaks and troughs in the 
ripple. 

• For each peak or trough determine the significance of the ex- 
cess or deficit in count rate in that region in the actual count rate R 
compared to the count rate that would result from the underlying 
surface brightness profile B. 

• Retain peaks or troughs only if their significance is greater 
than some threshold a; typically a = 3, is chosen. 

At the end of this process, we have a list of radii that mark the 
end-points of regions containing significant excesses or decrements 
in emission. However such regions do not necessarily constitute 
ripples. In order that we see a region as a ripple there must be some 
adjacent regions of positive or negative excess emission. The effect 
of varying the number of adjacent regions required to constitute a 
detection is examined in section [3".4. 1 1 

3.3 Generating model counts profiles 

In order to use the procedure above to investigate the dependence 
of detection time on the cluster and ripple properties, we require a 
method for generating realistic count rate profiles for clusters con- 
taining ripples. To generate these profiles, we start with functional 
forms for the underlying surface brightness profile, B, and a model 
for the modulus of the Fourier transform of the fractional residuals 
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Figure 4. Positive frequency half of the Fourier transform of the fractional 
ripples in the Perseus cluster (points) and in our model (thick red line). 

(i.e. the square root of the power spectrum), |J?(iJ)|. To generate 
the spatial surface-brightness ripples, we assign a random phase to 
each frequency component in the spectrum of the residuals, enforc- 
ing the condition that JP(R)(k) = ^'(i?) t (-fc), so that R(r) is 
real. The resulting surface brightness profile S = B(l + R) gives 
a random wave with the desired power spectrum. We initially gen- 
erate such a surface brightness profile with a high spatial resolution 
and then average over the radial bins appropriate to an observation 
to produce the observed counts profile. Poisson noise appropriate to 
the observation time is added when this surface brightness profile 
is converted to a count rate. By using different random frequency 
components, we are able to generate multiple spatial profiles and 
so average the estimated detection times over multiple clusters with 
the same ripple power spectrum. 

For our subsequent analysis, we often wish to compare re- 
sults from model clusters to those obtained for the Perseus clus- 
ter data. Therefore we construct a fiducial model for the ripple 
power spectrum which is based on the spectrum observed in the 
Pers eus cluster. We use the surface-brightness data shown in Fig. 
3 of lSanders & Fa bian (2007). This profile was generated in a sec- 
tor with an opening angle of approximately 0.13 x 2n, which we 
assume throughout this paper is a typical angle over which ripples 
will be coherent enough to allow a profile to be constructed from 
circular annuli. The model spectrum along with the actual Perseus 
spectrum is shown in Fig. [4] The model is of the form: 

\-9{R{r))\ = (1 - e~ ah ) * (L (fc) + Li(k) + B * e~ ck ) (12) 

where k is defined as 1/A and Lo and L\ are Lorentzian functions 
representing the main peaks in the surface brightness profile: 

AT 

L{k) = 2 7 r((z-* ) 2 +rV4) (B) 

Owing to the large number of free parameters in the model, 
the fit to the Perseus model was performed by eye, aiming to recre- 
ate the main features of the profile rather than provide an exact 
representation. 
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Figure 5. Number of ripples against observation time for the Perseus cluster 
with a detection requiring one (solid line), two (dashed line), three (dotted 
line) and four (dash-dotted) line features significant at the 3<r level. 

3.4 Results 

3.4.1 Observation Time to Detect Ripples in Perseus 

As a test of our method, we calculate the time required to detect 
ripples in the Perseus cluster. The projected beta models were fitted 
over the radius range 10 — 80kpc, and only ripples in this range 
were counted. 

Fig- H] shows the expected number of ripples detected against 
observing time for the Perseus cluster data. Each line represents a 
different choice for the number of significant excesses/decrements 
in emission needed to count a feature as a ripple. To account for 
the random nature of the Poisson noise added to our degraded sur- 
face brightness profiles, we have averaged the results over several 
realisations of the counts profile. For the case where 2-3 adjacent 
features are needed to constitute a ripple, we start to see several 
ripple-like features detected to 3u significance in about 40— 100 ks, 
broadly compatible with the bounds 20 ks < t < 200 ks imposed 
by the observed (non)detection times for ripples in Perseus. 

Fig. [6] shows the average number of ripples detected against 
time for a set of model clusters with the same ripple power- 
spectrum as Perseus. Again the different lines represent a different 
number of adjacent features needed to constitute a ripple. The num- 
ber of ripples detected as a function of time is comparable between 
the Perseus data and the model. In the case where a single feature 
counts as a several ripples are detected in just over 20 ks, but in 
the more realistic case where 2 — 3 adjacent features are needed to 
identify a ripple detections take 40 — 100 ks. 

3.4.2 Variation of Detection Time with Total Flux 

In order to determine how the time to detect a ripple varies with 
various parameters of the underlying surface brightness and the rip- 
ple structure, we use a binary search over the observation time. At 
each observation time, several realisations of the model cluster are 
constructed and we use the algorithm discussed above to search 
for ripples in the region r > lOkpc. At each time, we determine 



Figure 6. Number of ripples detected against observation time for an av- 
erage over model clusters with ripples having the same power spectrum as 
the Perseus cluster, with a detection requiring one (solid line), two (dashed 
line), three(dotted line) and four (dash-dotted) line features significant at 
the 3cr level. 



whether the median number of ripples detected is higher or lower 
than some threshold. The time for the next step is chosen as half 
time between the shortest time in which a detection has been made 
and the longest time in which no detection has been made. Initially 
we assume that ripples will always be detected in 10 7 s and will 
never be detected in 10 4 s (violation of these assumptions will lead 
to points lying just above the minimum or just below the maximum 
time, which will be obvious in the results). We typically consider 
an average of three ripple to be a detection and require at least two 
adjacent significant features to define a ripple. 

We use the same Perseus-like model cluster as before, with the 
same radial bins. This is necessarily an approximation; in reality 
one might choose larger bins for a fainter cluster. At each observa- 
tion time, 10 model clusters are generated to calculate the average 
number of ripples. 

Fig.|7]shows two independent determinations of this variation 
in the required detection time with the total flux of the cluster. The 
error bars represent the accuracy at which the binary search is ter- 
minated, that is, the difference between the maximum observation 
time for a non-detection and the minimum for a detection. To gauge 
the uncertainty in the detection time predicted by this method, we 
may simply compare points from the independent runs; these seem 
to agree well, with the discrepancies being similar in magnitude to 
the plotted error bars, indicating that we are terminating the binary 
search at a sensible accuracy. 

The detection time calculated for a Perseus-like cluster, with 
a flux of 6.75 counts per second inside 100 kpc is ~ 40 ks. 



3.4.3 Variation in the Detection Time With Ripple Amplitude 

Fig. [U shows to independent calculations of the detection time for 
various ripple amplitudes. The scaling of the detection time with 
ripple amplitude is roughly t ~ amplitude -2 , the same as lowest 
order in / of equation lllOt . 
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Figure 7. Two runs of the time required to detect Perseus-like ripples 
against the cluster flux in the detection region inside 100 kpc. The dashed 
line shows t ~ flux -1 , as expected from equation 4 10t . and is not a fit but 
is normalised for a detection time of 40 ks for the Perseus cluster (with a 
flux of 6.75 ct s _1 inside 100 kpc in the appropriate sector). The error bars 
account only for the accuracy of the binary search, not for the uncertainty 
in the time needed to make the detection. 
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Figure 8. Two runs of the time required to detect Perseus-like ripples 
against the ripple amplitude. The dashed line shows t ~ amplitude -2 
scaled so the detection time at the Perseus ripple amplitude is 40 ks. 



Figure 9. Two runs of the time required to detect Perseus-like ripples 
against the ripple wavelength. The wavelength shown corresponds to the 
highest frequency Lorentzian peak in the model. 



also affect the other frequency components. Also, as discussed in 
section [2] a ripple of smaller wavelength and equal amplitude in 
projection corresponds to a larger amplitude ripple in 3D. 

Fig. [9] shows the variation in detection time with the ripple 
wavelength. The detection time appears to be almost independent 
of wavelength above wavelengths of about 8 kpc. This result is in- 
consistent with a naive prediction; the number of counts in the peak 
of a ripple is N ~ SA where S is the surface brightness at the rip- 
ple radius and A is the area covered by the ripple. Typically, we 
have S ~ r -1 so N ~ r -1 rA ~ A, and so we expect the time to 
detect the ripple scales like ~ 1/A. To understand this discrepancy, 
we have run simulations using a simple monochromatic ripple pro- 
file, shown in Fig. [10] With monochromatic waves, the detection 
time more closely follows t ~ 1/A out to large ripple wavelengths. 
There is also some evidence that the beta-model underlying profile 
has larger detection times at radii above r ~ r colc than the scale- 
free power-law model. 



3.4.5 Variation in Detection Time with Underlying Cluster 
Properties 

Figure [TTI shows the variation in detection time with the (3 param- 
eter of the underlying beta model. For reasonable values of (3 there 
is little evidence of a substantial change in the detection time, al- 
though the time becomes large for very large (3 values. Variations 
in the core radius of the beta model similarly make little difference 
to the results. 



3.4.4 Variation in Detection Time with Ripple Wavelength 

In order to determine the variation in detection time with the rip- 
ple amplitude, we took the fiducial model and altered the values 
of xo in the two Lorentzian peaks, keeping the ratio of the two 
wavelengths constant, whilst leaving the rest of the spectrum unal- 
tered. This is inevitably an over-simplification; the underlying rip- 
ple wavelength is determined by the inflation process, which will 



3.4.6 Variation in Detection Time With Redshift 

Fig. [12] shows the variation of detection time with redshift at con- 
stant flux. This means that redshift variations correspond to a 
rescaling of the cluster's physical dimensions. For the cluster pa- 
rameters and range of redshifts shown, the detection time does not 
vary significantly, however there will be a point where the wave- 
length of the ripple is comparable to the size of the angular bins at 
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Figure 10. Time to detect monochromatic ripples against the ripple wave- 
length for underlying surface brightness profiles of beta model form with 
fcore = 10 kpc (circles) and with a power-law profile (green stars). The 
profiles are scaled so that the detection times at r = 10 kpc are similar. 
The blue dashed line indicates t ~ A -1 . 
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Figure 11. Mean time for detection against underlying cluster beta param- 
eter. 



the cluster redshift. After this point, which will occur at smaller red- 
shift for smaller wavelength, the ripples will become significantly 
more difficult to detect. 
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Figure 12. Mean time for detection against cluster redshift, assuming the 
cluster flux is constant, independent of redshift. 



4 DISCUSSION 

4.1 Time Needed to Detect Ripples 

Assuming that the different parameters in section 13.41 may be 
treated independently, we may derive a simple analytical approx- 
imation for the time to detect ripples in a given cluster. If we as- 
sume that the ripples are long enough wavelength that the scaling 
of detection time with ripple wavelength may be ignored, the only 
strong dependencies seen are with flux and amplitude, which scale 
in a manner consistent with the simple analytical model; t ~ flux" 1 
and t ~ amplitude" 2 . Taking the normalisation from the Perseus 
model: 



t = 40 



/ 



/perseus 



ks 



(14) 



for a cluster of flux / and ripple amplitude a. Except for the wave- 
length scaling, this result is very close to the analytic prediction of 
dlOl ), with the normalization increased by a factor ~ 3.5. 



4.2 Detectability of sound waves in nearby clusters 

To investigate the feasibility of detecting ripples in clusters other 
than Perseus using the current generation of X-ray satellites, we 
have used equation d!4t to determine the regions of a/ap axui , 
///perseus parameter space accessible to Chandra observations of 
different lengths. We have then estimated appropriate values for 
a/aperseus and ///perseus for several nearby clusters. 

To estimate ///perseus, we need the flux in ct s" 1 over the 
region where ripples could be detected in the cluster. We assume 

this corresponds to the region 2?"biibble centre Sbubble centre ? 

where 

Hubble centre is the a verage radius of the bub ble fr om the centre of th e 
cluster taken from lDunn & Fabianl J2004l) and lDunn et"all d2005l) . 
Making different assumptions about where the flux should be mea- 
sured does not change our results substantially. To measure the flux, 
we use raw Chandra events files taken from observations in the 
Chandra archive. In all cases except Virgo, the observations are 
taken with the ACIS-S instrument; for Virgo we scale from ACIS-I 
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to ASIS-S by increasing the count rate by a factor 1.5 appropriate 
for a 3keV plasma. The energy range is limited to 0.5 — 7keV. 
In each case, obvious point sources in the region of interest were 
removed by eye and a blank-sky background scaled to the correct 
exposure time and area was subtracted. 

To estimate the amplitude of the ripples, we use the relation- 
ship between the power in a sound wave and its pressure amplitude 
lLandau & LifshitzUl959l) : 



-Pwave = 47T7-' 



M 

pc 



(15) 



Assuming that the effect of projection is to reduce the projection 
reduces the surface brightness perturbation so 8S/S = nSn/n- as 
discussed in Section[2]- the amplitude of surface brightness fluctu- 
ations is related to the power in the wave as: 



r5S 
S 



kP 2 



(16) 



To estimate the wave power in this expression, we make 
use of two approaches; an estimate that P W ave = Jcavitv based 
on the cavity power using th e buoyancy timescale iBirzan et al.l 
(2004); Dunn & FabiarJ {2004); Raff ertv et all §006) and an esti- 
mate Pwave = icooi ( 1 - 2r bu bbie / ^cooi ) • The factor ( 1 - 2r coo i / r bu bbie ) 
is to account for the fact that the very central region of the cluster is 
not likely to be heat ed by sound waves but instea d by weak shocks 
and cavity heating I McNamara & Nul sen 2007). To estimate the 
projection factor k, we assume that the dominant ripple wavelength 
is ~ 2rb i.e. twice the average bubble dimension; this is a good ap- 
proximation for the Perseus cluster, k is then calculated using the 
procedure in section[2] We take r to be 5rbubbie centre- 
Using this model, wi t h the cluste r and bubble properties pre- 
sented in founn & Fabianl [2004, 2006) assuming that the cavity 
power is 4pV appropriate if the bubbles are filled with a fully- 
relativistic plasma, we predict the value of 8S/S for several local 
clusters, which we compare to SS/S = 0.05 for the Perseus clus- 
ter. Fig.fOJshows the position of these clusters on a flux-amplitude 
plane, including values for the Perseus cluster both calculated from 
this model and from the actual observations. The corresponding de- 
tection times with Chandra are indicated by the shaded regions. 

The above analysis suggests that detecting ripples in several 
local clusters may be possible with the Chandra satellite. In par- 
ticular Centaurus, Abell 1795 and Abell 2199 are promising candi- 
dates for detections in around 1 Ms of observation time. However 
it is important to note that there is considerable uncertainty in the 
wave amplitude determined for each cluster, with the two calcu- 
lated estimates giving up to an order of magnitude difference in the 
estimated amplitude. In Perseus it is apparent that the cooling lu- 
minosity is a better estimator of the wave amplitude than the cavity 
power, which is likely true for other clusters if the power in sound 
wave is closer to the time-average heating power than the instanta- 
neous cavity power. 

In addition to the power estimates, there are also uncertainties 
associated with the ripple wavelength and the radius at which it is 
assumed that the detection will be made. For a cluster like Centau- 
rus which is currently under-heated by the cavities, the final size of 
the bu bbles is likely greater than the current size; iDunn & Fabianl 
(2006) estimate rb/Vb, max = 0.75 for Centaurus, although it is no- 
table that many clusters in their sample have rb/rb,max > 1. indi- 
cating the cavities grow larger than needed to offset the cooling. 
Increasing the bubble radius will increase the wavelength of the 
ripples generated by the cavity inflation, reducing the projection 
factor for the ripples. For Centaurus, increasing the wavelength by 
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Figure 13. Flux and predicted ripple amplitude of some local clusters com- 
pared to Perseus. Blue plusses are based on the X-ray luminosity within 
the cluster cooling radius, red crosses are an estimate based on the cavity 
power. The black circle indicates the observed position of the Perseus clus- 
ter. Shaded regions indicate the predicted detection times for sound waves 
in the clusters using Chandra. 



a factor 1.5 would reduce the time needed to detect the waves by 
a factor of ~ 2.25 to a few hundred kilo-seconds. Conversely, the 
most prominent bubbles in Abell 2199 are likely detached from the 
central source and buoyantly rising in the ICM. Therefore their di- 
mensions likely overestimate the ripple wavelength and excluding 
the region inside 2rbubbic centre probably underestimates the cooling 
luminosity. This may explain the large discrepancy between the 
cooling luminosity and cavity luminosity based estimates of the 
amplitude in this system. 

Detection times may also increase if the cluster is in a regime 
where the wavelength of the ripples is important; for example if it 
has astrongly monochromatic ripples so the detection time scales 
as~ 1/A. 

Given the large uncertainty on the amplitude of the waves, 
and the strong dependence of the detection time on the amplitude, 
it may be possible to see ripples in some local clusters with cur- 
rently available data. In Abell 2199, an isothermal shock has been 
seen in a 33 ks exposure, but an unsharp-mask analysis revealed no 
further ripple-like features. For Centaurus, 200 ks of Chandra data 
are available, and have been analysed for ripples by Sanders et al. 
(in prep.). 

It is also important to recognise the limitations of our predic- 
tions compared to the procedure that has been used to find ripples 
in practice. The starting place to locate ripples in practice will be an 
image of the cluster such as Fig.Q]in which the small-scale features 
have been brought out using a technique such as unsharp-masking 
or a Fourier-filter. Assuming ripples exist in the cluster, whether 
they are seen in the resulting image will depend substantially on 
the degree of coherence of the ripples and the complexity of the 
radial power spectrum of the ripples. Ripples with a high degree of 
coherence might be discovered in much less time than suggested by 
our analysis whilst those with a limited degree of coherence might 
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not be noticed even in very long exposures. Since coherent features 
are likely easier to identify with better spatial resolution this sug- 
gests ripples will be easier to detect in low-redshift clusters than in 
higher redshift clusters even where their observed flux is similar. 

Related to this issue is that of the threshold number of standard 
deviations from the background count rate needed to class a feature 
as a ripple. In our analysis, we have assumed each feature needs to 
be detected to 3a. In practice, detecting a large number of clearly 
ripple-like features to less than 3a may be a more convincing de- 
tection than a detection in which a few blobs are detected to higher 
significant, but the intermediate structure is not clearly ripple-like. 

Given these limitations, extreme caution must be made in pre- 
dicting exact detection times for ripples and for interpreting non- 
detections as indicating that no ripple-like structures exist in a given 
cluster. 



5 CONCLUSION 

In order to understand the role played by sound waves in distribut- 
ing energy in cluster cores, it is essential to study these waves in a 
variety of cluster systems. To understand the requirements for such 
a study, we have investigated the effect of projection on reducing 
the observed wave amplitude compared to the intrinsic amplitude 
in a range of cluster atmospheres, and we have calculated the de- 
tection times for waves in a number of cluster environments. 

Projection of ripples in the emissivity profile substantially re- 
duces the amplitude of the surface brightness ripples. The magni- 
tude of this effect is a strong function of the ripple wavelength, but 
depends little on the underlying atmosphere properties. If the wave- 
length of the ripples is correlated with the dimensions of the cavity, 
this implies systems with larger cavities may be more promising 
targets for the detection of ripples. 

By constructing an algorithm for detecting ripples in model 
data, we have shown that the detection time for Perseus-like sur- 
face brightness ripples is critically determined by two features of 
the cluster - the total flux and the ripple amplitude. Other factors 
such as the underlying cluster properties have a much more limited 
effect. Applying our results to nearby clusters, and assuming that 
the ripple amplitude is sufficient to heat the cluster, we estimated 
detection times for ripples using the Chandra satellite. 

These detection times suggest that a selection of nearby, bright 
clusters may contain ripples detectable in around ~ 1 Ms of Chan- 
dra observation time, although there is considerable uncertainty 
brought about by uncertainties in the expected amplitude of the 
ripples. In cooler clusters such as Virgo that require less heating, 
and so are expected to have smaller ripples, the detection times are 
likely prohibitively long with Chandra but should be well within 
reach of XEUS which promises one to two orders of magnitude 
greater effective area. 
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